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ABSTRACT 

Multiplex analytical systems that allow detection of 
multiple nucleic acid targets in one assay can provide 
rapid characterization of a sample while still saving 
cost and resources. However, few systems have 
proven to offer a solution for mid-plex (e.g. 10- to 
50-plex) analysis that is high throughput and cost ef- 
fective. Here we describe the combined use of fluor- 
escence color and melting temperature (T m ) as a 
virtual 2D label that enables homogenous detection 
of one order of magnitude more targets than current 
strategies on real-time polymerase chain reaction 
platform. The target was first hybridized with a pair 
of ligation oligonucleotides, one of which harbored 
an artificial sequence that had a unique T m when 
hybridized with a reporter fluorogenic probe. The 
ligated products were then amplified by a universal 
primer pair and denatured by a melting curve 
analysis procedure. The targets were identified by 
their respective T m values in the corresponding fluor- 
escence detection channels. The proof-of-principle 
of this approach was validated by genotyping 15 
high-risk human papillomaviruses and 48 human 
single-nucleotide polymorphisms. The robustness 
of this method was demonstrated by analyzing a 
large number of clinical samples in both cases. The 
combined merits of multiplexity, flexibility and simpli- 
city should make this approach suitable for a variety 
of applications. 

INTRODUCTION 

Multiplex nucleic acid detection is a preferred approach in 
many applications such as screening for risk factors, 



typing of pathogen and detection of genetic disease pre- 
disposition. Multiplex essentially means the ability to sim- 
ultaneously ask more than one question about a sample at 
the same time, or the ability to conduct more than one test 
on a sample in the same reaction. Therefore, multiplex 
assays not only facilitate timely and cost-effective detec- 
tion, but they also offer more information in a sample 
than singleplex detections (1). Various multiplexing tech- 
niques have been described using measurement parameters 
detected on certain platforms. For example, multiplex 
polymerase chain reaction (PCR) exploits the length dif- 
ference of amplicons to identify targets through gel elec- 
trophoresis. Similar strategy can be found in multiplex 
ligase-dependent probe amplification (MLPA), where up 
to 40 nucleic acid sequences could be detected in a single 
capillary electrophoresis (2). Array systems, such as planar 
array and a variety of particle arrays, have been developed 
for suggested use in high or medium density multiplexed 
assay (3,4). In addition, molecular weight has been used 
as tags for single-nucleotide polymorphism (SNP) gen- 
otyping or mutation detection on the mass spectroscopy 
platform (5,6). So far, few DNA testing platforms have 
proven to offer a solution for mid-plex analysis that is 
high throughput, ease-of-use and cost effective. 

The emergence of real-time PCR has revolutionized 
nucleic acid detection in many respects owing to the 
reduced detection time and low risk of amplicon contam- 
ination. As a homogeneous detection by nature, real-time 
PCR distinguishes itself from the above-mentioned 
approaches in that the amplified products are detected in 
a closed-tube without post-PCR manipulations (7). 
Multiplex real-time PCR can be performed using a 
real-time thermocycler that has more than one detection 
channel, and thus allows different fluorophore-labeled 
probes to be simultaneously detected in one reaction. 
Because the number of detection channel, fmorophore 
and target are equivalent to each other in a classical 
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detection strategy, a real-time thermocycler equipped with 
four detection channels could accommodate a quadruplex 
detection. Higher multiplexed assays on a real-time PCR 
can be achieved using melting temperature (T m ) as a 
second dimension of multiplexing (8). For example, a 
recent study demonstrated that using four differently 
labeled probes, 27 different species of mycobacteria 
could be identified in one reaction (9). However, this 
T m -based multiplexity was restricted to targets differed 
in one variant region encompassed by a single pair of 
primers. Alternatively, a multicolor combinatorial probe- 
coding method has been proposed to detect multiple 
targets in a single reaction (10). The limitation of this 
strategy is that it is difficult, if possible, to detect 
multiple co-existing targets in one reaction. 

In the present study, we sought to develop a new multi- 
plex approach to overcome some of the aforementioned 
difficulties while maintaining the homogeneous nature of 
real-time PCR. This approach is achieved by the combina- 
torial use of the fluorophore and T m as 2D label, which 
enables an order of magnitude increase in multiplexity on 
the real-time PCR platform. As a proof of principle, the 
2D label was first used to establish a 16-plex assay to 
genotype 15 high-risk human papillomaviruses (HPVs). 
This assay aimed to test whether the 2D label could be 
used to detect one or several among all possible targets in 
a single reaction. The robustness of this assay was 
validated by a comparison study with 517 clinical 
samples. In a second assay, the 2D label was used to es- 
tablish a 96-plex assay for genotyping of 48 forensic SNPs. 
The purpose of this assay was to determine whether the 
2D label could be used to detect all possible targets in a 
single reaction. The established assay was validated with 
100 human DNA samples in a head-to-head comparison 
with a second commercially available method. Moreover, 
the utility of this assay was demonstrated in assessment of 
hematopoietic chimerism after bone marrow transplant- 
ation. The result from both assays demonstrated that the 
2D label can be used as a universal label to detect the 
existence of one or a few among all possible targets or 
even simultaneous existence of all possible targets in one 
reaction. 



MATERIALS AND METHODS 

Construction of a library containing 50 2D labels 

A library containing 50 2D labels was constructed by 
using five fluorophores (X = 5) and 10 T m tags (Y = 10) 
(Supplementary Table SI). The five types of fluorophores 
were FAM, HEX, ROX, CAL Fluor Red 635 and Quasar 
705. The sequences of the fluorophore-labeled probe were 
artificially generated and are not homologous to any 
known species. For each fluorophore-labeled probe, 10 
T m tags with predicated T m values ranging from 40°C to 
80°C with 1°C-5°C intervals were designed to hybridize 
with it. The T m tag of the highest T m is fully complemen- 
tary to the probe. Other tags consist of different number 
of mismatched bases that are introduced to confer a 
gradual decrease in T,„ values. All synthetic oligonucleo- 
tides were prepared at the polyacrylamide gel 



electrophoresis purity by Sangon, Inc. (Shanghai, China) 
and the fluorophore-labeled probes were from Biosearch 
Technologies, Inc (Novato, CA, USA). 

The predicated T m values were obtained by T,„ Utility 
vl.3 (BioFire Diagnostics Inc., Salt Lake City, UT, USA) 
using the following conditions: 200 nM probe, 200 nM 
target, 3mM MgCl 2 and 200 uM deoxyribonucleoside 
triphosphates (dNTP), where the probe was the 
fluorogenic probe and the target was the T m tag. The 
actual T m values were measured by melting curve 
analysis on a Rotor-Gene™ 6000 real-time analyzer 
(Corbett Research, Mortlake, Australia) in a 25-ul 
solution of 10 mM Tris-HCl (pH 8.6) containing 50 mM 
KC1, 3mM MgCl 2 , 200 uM dNTP, 200 nM fluorogenic 
probe and 200 nM T m tag. The melting curve analysis 
(MCA) procedure started from 95°C for 1 min, 35°C for 
2min followed by raising the temperature from 40° C to 
85°C at l°C/step with a 5s stop between each step. 
Fluorescence intensity was recorded at each step in each 
of the five detection channels. The data obtained were 
plotted as fluorescence intensity versus temperature as 
well as the negative derivative of fluorescence intensity 
with respect to temperature. T„, values were identified 
from the peak position of the melting curve. 

Typing of 15 high-risk human papillomavirus genotypes 
Step 1: Hybridization and ligation 

The reaction was performed in a 10-ul solution containing 
20 mM Tris-HCl (pH 7.6), 25 mM sodium acetate, 10 mM 
magnesium acetate, 1 mM dithiothreitol, 1 mM NAD + , 
0.1% Triton X-100 containing 30 type-specific ligation 
oligonucleotides (500 pM each, Supplementary Table 
S2), two ligation oligonucleotides for the internal control 
(IC) (250 pM), 1 unit of Taq DNA ligase (New England 
Biolabs, Bejing, China), 5 x 10 2 copies SUC2 gene (IC) 
and 5ul of either artificial plasmid DNA templates ex- 
tracted from E. coli DH5oc by a plasmid mini kit 
(Omega Bio-Tek, Inc. Norcross, GA, USA) or DNA ex- 
tracted from clinical specimens by a kit supplied with the 
HPV Geno Array test kit (Hybribio Ltd., Chaozhou, 
China). Clinical specimens were provided as on-shelf 
coded samples that had been collected previously for diag- 
nosis purposes. The reaction was performed in a T3 
thermocycler (Biometra, Gottingen, Germany) at the fol- 
lowing temperatures: denaturation at 95°C for 5 min; 10 
cycles of 95°C for 1 min, 70°C for 1 min, 68°C for 1 min, 
66°C for 1 min and 64°C for 3 min. 

Step 2: PCR/MCA 

The reaction was performed in a 25-ul solution containing 
10 mM Tris-HCl (pH 8.6), 50 mM KC1, 2.5 mM MgCl 2 , 
5% (V/V) glycerin, 1 U Taq DNA polymerase, 1 pmol 
limiting primer, 20 pmol excess primer, 5 pmol of 
fluorogenic probes for 2D label (Supplementary Table 
S2) and 2ul reaction product from step 1. PCR/MCA 
was performed in the Rotor-Gene™ 6000 at the following 
temperatures: denaturation at 95°C for 3 min, 42 cycles of 
95°C for 10 s, 56°C for 10 s and 72°C for 20 s, followed by 
a temperature increase from 40°C to 85°C, l°C/step, with 
a 5s stop between each step. Fluorescence intensity was 
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measured in four detection channels: FAM (510 nm), 
HEX (555 nm), ROX (610 nm) and CAL Fluor Red 635 
(660 nm). 

Genotyping of 48 SNPs 

Step 1: Hybridization and ligation 

The reaction was performed in a 10-ul solution containing 
20 mM Tris-HCl (pH 7.6), 25 mM potassium acetate, 
10 mM magnesium acetate, 1 mM dithiothreitol, 1 mM 
NAD + , 0.1% Triton X-100, 144 SNP-specific ligation 
oligonucleotides (800-8000 pM each), 1 unit of Taq 
DNA ligase and 5 ul of human genomic DNA extracted 
from anonymous blood samples obtained from Xiamen 
Blood Center by QIAamp DNA Mini Kit (Qiagen, 
Hilden, Germany). The reaction was performed in a T3 
thermocycler at the following temperatures: denaturation 
at 95°C for 5min; 10 cycles of 95°C for 2min, 68°C for 
2min, 66°C for 2min, 64°C for 2min, 62°C for 2 min and 
60°C for 2 min. 



The design of 2D label assay 

The design of 2D label assays is illustrated in Figure 2. 
Two ligation oligonucleotides are designed for one target. 
The left ligation oligonucleotide contains an upstream 22- 
nt universal priming sequence at the 5'-end, a 4-nt stuffer 
sequence, the T m tag sequence and a 26^10-nt target- 
specific sequence at the 3'-end. The right ligation oligo- 
nucleotide contains a 30-50-nt target-specific sequence at 
the 5' phosphorylated end and a downstream 21-nt uni- 
versal priming sequence at the 3'-end (Figure 2a). The T,„ 
tag in the left ligation oligonucleotide is fully or partially 
complementary to a particular fluorogenic probe and thus 
creates the specific X and Y values of a 2D label for the 
target. Melting temperatures of all probe target-specific 
sequences are 65°C-70°C under conditions of ligation. 
For the design of the hybridizing parts, sequence informa- 
tion available from the public databases (www.ncbi.nlm. 
nih.gov) was used. 



Step 2: PCR/MCA 

The reaction was performed in a 25-ul solution containing 
75 mM Tris-HCl (pH 9.0), 20 mM (NH 4 ) 9 S0 4 , 50 mM 
KC1, 3mM MgCl 2 , 0.01% (V/V) Tween 20, 1U hot 
start Taq DNA polymerase, 200 liM dNTP, 1 pmol 
forward primer, 20 pmol reverse primer, 5 pmol of 
fluorogenic probes for 2D labeling and 1 id of reaction 
product from step 1. PCR/MCA was performed in a 
Rotor-Gene™ 6000 at the following temperatures: de- 
naturation at 95°C for 3 min, 50 cycles of 95°C for 5 s, 
56°C for 10 s and 75°C for 15 s, followed by a temperature 
increase from 35°C to 90°C, l°C/step, with a 5 s stop 
between each step, and fluorescence intensity was 
measured in four detection channels, as described above. 



RESULTS 

The working principle of 2D label 

In 2D label, the type of different fluorophores (X) defines 
the first dimension and the value of resolvable T m (Y) 
defines the second dimension. The total number of labels 
is equal to the product of X multiplied by Y. Figure la 
shows an example that combines five fluorophores (X = 5) 
and six T,„s (Y = 6), resulting in 30 different labels. To 
create a 2D library, a series of oligonucleotides, called Tm 
tags, are used to form hybrids with a fluorophore-labeled 
probe, resulting in different characteristic T„,s from this 
fluorophore. As shown in Figure lb, six oligonucleotides 
and one FAM-labeled probe yield six 2D labels. 

In a spectrofluorometric thermal cycler, the number of 
differently colored fluorophores that can be used is 
determined by the number of detection channels, which 
usually ranges from two to six, depending on the instru- 
ment. On the other hand, the number of T m s that can be 
reliably distinguished from a melting curve varies from 6 
to 10, depending on the temperature resolution of the in- 
strument. Consequently, the maximal number of 2D labels 
can be detected is 12-60 with instruments currently 
available. 



Fluorophore 




2D Label (X=5) 



(b) 



Probe 5'-FAM-(GGTCCTTCATCGCTCAGCCTTCACCGG)-BHQ1-3' 2D Label 

(X=FAM, Y=T m ) 

Target 1 5'-(GGTCCTTCATCGCTCAGCCTTCACCGG)-3' (T m =77.5°C) (FAM, 77.5) 

Target 2 5'-(GGTCCTTCATCGCTCAGCCATCACCGG)-3' (T m =74.0"C) (FAM, 74.0) 

Target 3 5'-(GGTCCTTAATCGCTCAGCCATCACCGG)-3' (T„=67.8X) (FAM, 67.8) 

Target 4 5'-(GGTCCTTAATCGCTCAGGCTTCACCGG)-3' (T m =63.6"C) (FAM, 63.6) 

Target 5 5'-(GGTCCTTGATCGGTCAGCCATCACCGG)-3' (T m =58.0"C) (FAM, 58.0) 

Target 6 5'-( GGTCCTTGATGTCTCAACCTACACCGG)-3' (T„=49.5°C) (FAM, 49.5) 

Figure 1. The concept of 2D label and its composition, 
(a) Combination of five fluorophores and six T,„s to create 30 2D 
labels. The five fluorophores (shown in color as fluorescence emission 
wavelengths) are those detectable in a common real-time PCR machine. 
The six T,„s (shown as six melting peaks) are obtained from melting 
curve analysis in one detection channel, (b) The labeling fluorophore of 
the probe, FAM (X = 1) and the six T,„ tags (Y = 6) together form six 
2D labels. Among the six T„, tags, the one with the highest T,„ is fully 
complementary with the probe, all others possess base pair substitu- 
tions (underlined) that result in mismatched base pairs in the probe: tag 
hybrid. The T„, variations of the hybrid are determined by the number 
and position of mismatched base pairs present. 



e76 Nucleic Acids Research, 2013, Vol. 41, No. 7 Page 4 of 1 1 

( a ) P1 S T m tag H1 H2 P2 
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(1) Hybridization (2) Ligation 



/ ^ / 



(3) LATE PCR (4) Melting 




T m 



Figure 2. The design of 2D label assays, (a) Two unreacted ligation oligonucleotides in 2D label assays. The left ligation oligonucleotide consists of 
four segments: an upstream priming region (PI), a stuffer sequence (S), a T,„-tag sequence that serves as the identification label of the target and a 
hybridization region (HI) that is complementary to the target sequence in the genomic DNA. The right ligation oligonucleotide consists of two 
segments: a hybridization region (H2) that is complementary to the target sequence in the genomic DNA, and a downstream priming region (P2). (b) 
Flow chart of 2D label assays. (1) The hybridization regions of the two ligation oligonucleotides hybridize to the target sequence. (2) The two 
oligonucleotides are ligated to each other, producing a joined full-length sequence. (3) The ligated product is amplified with a universal pair of 
LATE-PCR primers. The LATE-PCR primer that is present in excess is shown as a thick arrow, and the LATE-PCR primer that is present in limited 
concentration is shown as a thin arrow. (4) Melting curve analysis in certain detection channel yields a T„, value identifying the target sequence. 



During the hybridization-ligation process, the two oligo- 
nucleotides hybridize to the genomic DNA sample. 
Because hybridization occurs before any amplification 
steps, no amplification bias can be introduced into the 
assay. Following hybridization, a DNA ligase seals the 
nick between the left and the right ligation oligonucleotides 
to form PCR templates (11) that can be amplified with a 
universal PCR primer pair designed for asymmetric Linear- 
after-the-exponential (LATE) PCR (12). The LATE PCR 
is performed in the presence of a fluorogenic probe using a 
real-time PCR machine. After thermal cycling, a melting 
curve analysis allows readout of the unique T,„ values dis- 
played by the hybrids formed between T m tags and their 
corresponding fluorescent probes in the respective 
fluorophore channels, and thereby producing the X and 
Y values for the targets that exist (Figure 2b). 

Multiple targets can be detected using different target- 
specific pairs of ligation oligonucleotides. These targets are 
differentiated by their T,„ tags. We have designed a library of 
50T m -tag sequences, 10 of which serve as cassettes that can 
be detected by one fluorogenic probe. These pre-optimized 
T„,-tag sequences can be incorporated into any ligation 
oligonucleotides for multiplex assays of target DNA se- 
quences of interest. These T„,-tag-probe combinations are 
compatible with all mainstream spectrofluorometric 
thermal cyclers (Supplementary Table SI). 

Detection of 15 high-risk HPV types 

We set out to validate the 2D label scheme by detecting part 
of all the possible targets in a multiplex assay, which was 



exampled by detecting 15 high-risk HPV types in one 
reaction. Figure 3a illustrates the flow chart of the detection 
procedure. In this reaction, 16 target-specific ligation oligo- 
nucleotide pairs were present, and each pair possessed a 
unique T,„-tag sequence. Four differently colored probes 
were present, the FAM-labeled probe detected amplicons 
generated from an IC target sequence; the three other 
amplicon probes each detected the amplicons generated 
from five different HPV target sequences. The identity of 
each target sequence was revealed by the combination of 
probe color and probe-amplicon T,„. The entire procedure 
comprised a 16-plex ligation, a singleplex four-color PCR 
and an MCA. Using artificial plasmid DNA as HPV tem- 
plates, the MCA results showed that each of the 15 HPV 
types was correctly identified by its 2D label (Figure 3b). 

In clinical samples, it is common for each HPV type to 
have various subtypes owing to the existence of poly- 
morphic nucleotides. This added complexity was ad- 
dressed by using a mixture of ligation oligonucleotides 
that cover all polymorphic sites, but harbor one 
common T„, tag, such that all the polymorphisms belong- 
ing to the same type were identified as a single HPV type 
(Figure 4). Co-infection by more than one HPV type also 
occurs clinically. Using dual infection of HPV- 16 and 
HPV- 18 as an example, we studied the ability to distin- 
guish the two types when they are simultaneously present 
in the same reaction. By varying their relative ratios, the 
results showed that as low as 1% of one type in the 
presence of 99% of the other type could be detected 
(Figure 5). 
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(a) 16-plex Ligation Singleplex PCR/MCA (&) 

(90 min) (100 min) 




Temperature (°C) 

Figure 3. The application of 2D label to the typing of 15 HPV types in a single reaction, (a) Flow chart of the assay. In the ligation step, 16 pairs of 
ligation oligonucleotides were used to detect 15 HPV genotypes and one IC. In the PCR/MCA step, the ligated products were amplified by one pair 
of primers in the presence of four differently labeled fiuorogenic probes. T,„ tags are shown in different colors (corresponding with the labeling 
fluorophores in the fiuorogenic probes), and the black bars within the T„, tags represent those substituted mismatched nucleotides, (b) Readout for 
HPV typing using the 2D label method. Each HPV type was identified by its unique 2D label, i.e. the combination of a fluorophore type and a T m 
value. To visually illustrate the T m values, the melting curves were normalized between 0 and 1, and the normalized value portions between 0.4 and 1 
are shown. FAM: green line; HEX: orange line, ROX: red line and C635: blue line. The numbers above each melting peak identify the 
HPV genotype. 



We then studied the concentration range of the 
assay using plasmid templates. The results showed that 
5 x 10 2 5 x 10 7 copies plasmid DNA per reaction could 
be repeatedly detected regardless of the genotypes 
(Supplementary Figure SI). 

The robustness of the 2D label-based assay was 
validated by comparing genotyping results of 517 
cervical swab samples with a commercial chip-based 
system (GenoArray, Hybribio, Chaozhou, China) in a 
blind manner. The overall agreement between our 
method and the GenoArray assay was 96% (495/517), 
while the type agreement between the two methods was 
greater than 99% (Table 1). These results demonstrated 
the suitability of our 2D label-based strategy for clinical 
applications. 

Genotyping of 48 forensic SNPs 

We further tested the 2D label strategy by establishing a 
second assay to genotype 48 forensic SNPs that are 
randomly distributed among 23 human chromosomes. 



Unlike the HPV typing assay in which only a portion of 
the 1 5 types were present in a sample, the genotyping of 48 
SNPs could have 48-96 alleles simultaneously present in 
one reaction. To set up this assay, we designed two left 
ligation oligonucleotides harboring different T,„ tags for 
each allele and one right ligation oligonucleotide for each 
locus. Thus, the ligation step of this assay was a 96-plex 
reaction containing 144 ligation oligonucleotides. To 
identify all the 96 alleles of 48 SNPs, 3 PCR reactions 
were used, each covering 32 alleles of 16 SNPs 
(Figure 6a). Thirty-two 2D labels were used for each 
reaction (Supplementary Table SI). After PCR/MCA, 
the melting curves (Figure 6b) obtained from the four 
channels were imported into the GeneMarker software 
(Softgenetics, State College, PA, USA) for data trans- 
formation. The function Melt Analysis, which provides 
data interpolation, baseline calibration, wavelet trans- 
formation, fast Fourier transformation and automatic 
peak detection with allele identification, was used to 
analyze the melting data and identify the allele of each 
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Standard ligation oligonucleotide Degenerate ligation oligonucleotide 




50 60 70 80 50 60 70 80 

Temperature (°C) 



Figure 4. The effect of using degenerate ligation oligonucleotides on the performance of the 2D label assay. Both a standard ligation oligonucleotide 
and a degenerate ligation oligonucleotide mixture were used to detect a standard HPV-16 target and a polymorphic HPV-16 target, respectively. The 
polymorphic site is at the penultimate nucleotide at the 3'-terminus of the left ligation oligonucleotide sequence. The template DNAs used were 
artificial plasmids with concentrations of 1 x 10 4 , 1 x 10 3 , 1 x 10 2 , 1 x 10 1 and 1 x 10° copies/ul, respectively. Note that the melting curve peak of the 
1 x 10" copies/ul (dashed lines) was lowest when a standard ligation oligonucleotide was used to detect the polymorphic HPV-16. Non-template 
controls (NTCs) are shown as grey lines. 



SNP. A typical result from Reaction 1 was shown in 
Figure 6b. 

We investigated the lowest amount of human genomic 
DNA required for these genotyping assays. The result 
showed that at least 20 ng DNA per assay was needed to 
achieve reproducible genotyping readout. Such starting 
genomic DNA amount required for our 2D label assays 
was similar to that required for MLPA (2). 

The 48-SNP genotyping assay was used to analyze 100 
human genomic DNA samples and the results were 
compared with a commercial singleplex real-time PCR- 
based genotyping assay (ForenSNPs, Zeesan Biotec, 
Xiamen, China). A 100% agreement was achieved 
between these two methods (Supplementary Table S3). 
As a clinical application, the 48-SNP genotyping assay 
was used to assess hematopoietic chimerism after bone 
marrow transplantation. A patient who had a sibling 
bone marrow donor possessed an identical human leuco- 
cyte antigen with the donor. However, SNP 14 was found 
to have a homozygotic disparity between the donor 
and the recipient and thus was chosen to track the level 
of donor granulocyte engraftment. The results showed 
that donor monocytes began to appear 10 days after 
transplantation and became dominant 20 days after trans- 
plantation (Figure 7). This result indicated that a semi- 
quantitative clinical analysis could be carried out by 
using 2D labels. 



DISCUSSION 

The 2D label is representative of two measurable param- 
eters, i.e. fluorescence and T„„ which can be detected in a 
homogeneous way on real-time PCR machines. It was 
realized by using artificial T m tags and corresponding 
fluorogenic probes. In this study, we used ligation 
reaction as the labeling procedure, which allowed the use 
of a single universal primer pair to amplify all ligated 
products. Such a labeling strategy enabled the transform- 
ation of the multiplex assays into a singleplex PCR detec- 
tion. The total number of the targets that could be 
detected in a 2D labeled assay is equal to the number of 
readily resolvable T m values multiplied by the number of 
fluorescence detection channels of a real-time PCR 
machine. For a standard 4-color real-time PCR machine 
that has a T m resolution of 5°C in the range of 45°C-85°C, 
the 2D label strategy could accommodate a 36-plex detec- 
tion. Higher plex levels could be achieved if the machine 
has more detection channels and/or a better temperature 
resolution. 

The multiplexing capability of 2D label-based assays 
fills the gap between multiplex real-time PCR and 
array-based assays. Such an intermediate complexity is 
especially useful in molecular assays of routine use. For 
example, diagnosis of sepsis (13), gastrointestinal infec- 
tions (14) and respiratory diseases (15) often involve 
dozens of pathogens of different types including viruses, 
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bacteria and parasites. Similarly, inherited diseases and 
tumors are often caused by a variety of mutations found 
in different exons or genes (16,17). To detect these 
multiple targets or variants, current real-time PCR-based 
assays have to use multiple reactions, and other methods 
with higher order multiplexity require complicated 
post-PCR manipulations or expensive detection 
platforms. 

Using this new 2D label concept, we have successfully 
established two multiplex assays and obtained comparable 
results with those well-established methods. The HPV 
typing assay involved simultaneous detection of 15 high- 
risk HPV types responsible for cervical cancer and one IC 
from a single sample in a single reaction. The method 
could detect intratype variations and multiple types in 
mixed infections, had a wide detection range and was 
highly concordant with a commercial array-based system 
when tested with 517 cervical swab samples. Most import- 
antly, this assay could analyze 96 samples in 190min 
without any post-PCR manipulations and the consumable 
cost per assay was close to a singleplex 4-color real-time 
PCR. By contrast, the comparison method could only pro- 
cess up to 30 samples in a working day (18). There are 
reported assays for 14 HPV types in a single reaction that 
could be finished within a shorter time. However, such 
assays could not give a full genotyping result (19,20). 
The SNPs genotyping assay was able to genotype 48 
forensic SNPs randomly distributed among the 23 pairs 
of human chromosomes from a single sample in a single 
reaction. This is a sharp contrast with a reported melting 
curve analysis method that could only genotype six 
forensic SNPs in one reaction on a six-channel real-time 
PCR machine (21). The accuracy of the assay was 
validated by analyzing 100 human genomic samples with 
a complete concordance with the comparison singleplex 
real-time PCR method. In addition, the robustness of 
this assay was demonstrated in the semi-quantitative de- 
tection of the level of donor granulocyte engraftment after 
bone marrow transplantation. 



Table 1. Comparison of the HPV genotyping results of 517 specimens between two assays 



HPV type 




Number of samples 


(2D label/GenoArray) 


Number of discrepant results 


Agreement (%) 


Kappa (95% CI) 




+/+ 


+/- 


-/+ 










16 


49 


2 


0 


466 


2 


99.6 


0.978 (0.962-0.994) 


IS 


8 


2 


0 


507 


2 


99.6 


0.887 (0.808-0.966) 


31 


17 


1 


0 


499 


1 


99.8 


0.970 (0.940-1.000) 


33 


12 


1 


0 


504 


1 


99.8 


0.959 (0.918-1.000) 


35 


3 


0 


0 


514 


0 


100 


1.000 (1.000-1.000) 


39 


11 


3 


0 


503 


3 


99.4 


0.877 (0.807-0.947) 


45 


2 


0 


2 


513 


2 


99.6 


0.665 (0.442-0.888) 


51 


5 


2 


2 


508 


4 


99.2 


0.710 (0.572-0.848) 


52 


37 


3 


2 


475 


5 


99.0 


0.931 (0.901-0.961) 


56 


6 


2 


0 


509 


2 


99.6 


0.855 (0.754-0.956) 


58 


34 


1 


0 


482 


1 


99.8 


0.984 (0.968-1.000) 


59 


8 


0 


0 


509 


0 


100 


1.000 (1.000-1.000) 


68 


14 


1 


0 


502 


1 


99.8 


0.965 (0.930-1.000) 


73 a 


1 


0 


0 


516 


0 


100 


1.000 (1.000-1.000) 


82" 


3 


0 


0 


514 


0 


100 


1.000 (1.000-1.000) 



"HPV-73 and HPV-82 were not included in the GenoArray technique. Their detection results obtained from the 2D assay were compared with the 

sequencing method. 

CI = confidence interval. 




Temperature (°C) 

Figure 5. Simultaneous detection of multiple HPV genotypes in mixed 
infections, (a) Two co-existing genotypes were detected in the same 
channel, (b) Two co-existing genotypes were detected in two different 
channels. Melting curves of artificial plasmid templates containing two 
HPV genotypes of various percentage ratios (0:100, 1:99, 3:97, 5:95, 
10:90, 20:80, 30:70, 40:60, 50:50, 60:40, 70:30, 80:20, 90:10, 95:5, 97:3, 
99:1 to 100:0) were tested. The overall template concentration was 10 6 
copies per reaction. NTCs are shown as black lines. 
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(a) 

96-plex ligation in one reaction (120 min) 



Singleplex PCR/MCA in 3 reactions (120 min) 



Reaction 1 



Reaction 2 



Reaction 3 



F1 FAM BHQ1 
SNP(2, 18, 34) — »-•- 

A 
T 




SNP (5, 21, 37) 
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■ C J 
■T 



R1 F2 FAM BHQ1 



R2 F3 FAM BHQ1 



R3 



HEX BHQ1 



HEX BHQ1 



HEX BHQ1 



^ SNP (9, 25, 41) 
SNP (11, 27, 43) 



SNP (10, 26, 42^ 

SNP (12, 28, 44) 
■G~ 



ROX BHQ2 



ROX BHQ2 



ROX BHQ2 



SNP (13, 29, 45) 
.A 

G 

SNP (15, 31, 47) 



SNP (14, 30, 46) 

.C 

G 

_SNP (16, 32, 48) 

■a y 



C635 BHQ2 



C635 BHQ2 



C635 BHQ2 



(b) _ 

SNPs No 1 



FAM 



HEX 



ROX 



C635 



2 3 4 
Alleles TCTATC GA 



5 6 7 8 
CTTC TC TC 



9 10 11 12 
C A G AG A G C 



13 14 15 16 
GAG CGCGA 




50 



60 



70 



80 



T" 

70 80 50 60 
Temperature (°C) 



Figure 6. The application of 2D label in genotyping of 48 SNPs. (a) Flow chart of the assay. Forty-eight groups of ligation oligonucleotides 
(containing two left ligation oligonucleotides and one right ligation oligonucleotide for each SNP) were used to target 96 alleles of the 48 SNPs. 
The allele types of SNPs 1-16 were given at the 3'-terminus of the left ligation oligonucleotides. T„, tags are shown in different colors (corresponding 
with the labeling fluorophores in the fluorogenic probes), and the black bars within the T,„ tags represent those substituted mismatched nucleotides. 
In the PCR/MCA step, the ligation products were amplified by three independent pairs of LATE-PCR primers (Fl/Rl, F2/R2 and F3/R3) in the 
presence of the four fluorogenic probes, (b) Typical readout of 16 SNPs in reaction 1 from one human genomic DNA sample. Upper panel shows the 
melting curves and the lower panel shows the exported results from GeneMarker software. The allele types of the 16 SNPs were identified according 
to their detection channels. SNPs 1-4: FAM; SNPs 5-8: HEX; SNPs 9-12: ROX; SNPs 13-16: C635. 
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Figure 7. SNP genotyping for assessment of hematopoietic chimerism after bone marrow transplantation. The leftmost panel shows that SNP 14 has 
a homozygotic disparity between the donor (blue line) and the recipient (red line). The other three panels display the time-dependent change from 
allele G to allele C of SNP 14 in the recipient. 



Melting temperature was proposed as a second layer of 
multiplexing more than a decade ago to serve to augment 
multiplexing capability (8). However, high order multi- 
plexing remains to be achieved by existing strategy in 
which the T m is obtained from the probe-target hybrid. 
This strategy has been widely used for detection of muta- 
tions in close proximity or in one amplicon (22-25). 
However, if the variants are in different amplicons, such 
a strategy will encounter difficulty in dealing with the 
multiple primer pairs and adjusting the T,„ values to 
have balanced amplifications and discriminable T„, 
values. Another problem comes from the unexpected mu- 
tations (26) and polymorphism (27) that might occur in 
the probe-binding region of the targets. Such mutations 
and/or polymorphisms might give similar T„, value with 
the variant to be detected and cause false positive results. 
One mutation/polymorphic site may be easy to deal with 
but if multiple mutations/polymorphic sites are present, 
discrimination of variants by one probe becomes a 
challenging task. By comparison, the 2D label strategy 
uses a multiplex ligation reaction followed by a universal 
PCR, a well-established procedure that can detect dozens 
(2) and thousands (28) of loci from different amplicons. 
The presence of polymorphism in the target would exert 
no impact on the T,„ predefined for the target because the 
T,„ value was obtained from the T„, tag-probe hybrid 
rather than from the target-probe hybrid. Problems 
might be encountered if the polymorphic nucleotides 
occur and prevent ligation reaction. Such problems can 
be resolved by using degenerate nucleotide or universal 
binding base at the polymorphic sites or, if possible, 
changing the ligation region. Moreover, the use of T m 
tag could obviate the possible abnormal melting curve 
caused by secondary structure in the target, which other- 
wise might cause unexpected result in probe-target 
hybrid-based melting curve analysis. 

Compared with other multiplexing techniques, the 2D 
label approach has advantages in many respects. First, as 
a homogeneous assay by nature, its amplification and de- 
tection are carried out within a consecutive procedure and 
in a closed-tube format. Like in real-time PCR, this 
feature significantly simplifies the manipulations, reduces 
the risk for PCR product contamination and increases the 



assay throughput. Second, by using a common 2D label 
library, the assay is easy to be standardized in set-up, de- 
tection and readout. Such standardization helps save the 
overall, and especially labor, cost. In the two exampled 
assays, the amount of ligation oligonucleotides used in 
one reaction was 5fmol, and one standard commercial 
preparation for the ligation oligonucleotide (lOOnmol) 
can therefore be used for 20 million reactions! Therefore, 
once the assay was established, the ligation oligonucleo- 
tides cost was negligible. The only costly reagents are 
fluorogenic probes. However these universal probes can 
be prepared in bulk because the same fluorogenic probes 
can be used for different assays (e.g. for both the HPV and 
SNP assays). In this study, one fluorogenic probe was 
designed to detect all the T„, tags in one channel. It 
might be possible to use more than one fluorogenic 
probe in some applications, but the number of fluorogenic 
probes is still drastically smaller than the targets to be 
detected. Consequently, the material cost of the 2D 
label-based assay for a double-digit number of targets is 
close to multiplex real-time PCR assays for a single-digit 
number of targets. Finally, the 2D label-based assay can 
be performed with almost any real-time PCR machines, 
making it an open system readily adoptable for different 
applications. 

The proposed 2D label strategy for multiplex detection 
has several limitations worth further improvement. First, 
there is an extra ligation step before PCR in the 2D label 
assays when compared with the real-time PCR detection. 
This extra pre-PCR step increased the manipulations and 
lowered the throughput; however, it causes no risk of PCR 
amplicon contamination frequently encountered in 
strategies involving post-PCR manipulations. It is 
possible for this ligation reaction to be combined with 
PCR as reported previously (29), and therefore the add- 
itional step might be eliminated. Second, non-specific 
amplification might derive from the non-ligated oligo- 
nucleotides, producing false positive signals and lowering 
the overall analytical sensitivity. Although the non-specific 
signal can be largely reduced or eliminated by changing 
the ligation oligonucleotide design, one useful way is to 
include an IC in the system as exampled in the HPV assay. 
The amplification of the IC consumed the excess primer 
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and reduced opportunity for non-specific amplification in 
the absence of target. Another way to eliminate non- 
specific amplification is to add certain exonucleases to 
digest unligated oligonucleotides before amplification 
(30). It should be noted, however, that the concentration 
of the ligation oligonucleotides is extremely low compared 
with the universal primers, and the non-specific signals 
could be largely eliminated through design of the 
ligation oligonucleotide without digestion. Third, as an 
end-point detection, the 2D label-based assays are essen- 
tially non-quantitative and therefore restricting their use 
in screening rather than quantitative detection. However, 
as shown in the example of monitoring donor monocytes 
after bone marrow transplantation, a semi-quantitative 
detection could be achieved. Further research is now 
being undertaken to assess whether a combination of 
real-time PCR detection and melting curve analysis can 
improve the quantification ability. Finally, the analytical 
sensitivity of this method was relatively low as observed in 
the SNP genotyping assay, which might preclude its use in 
some forensic cases when extremely low amount of residue 
DNA was available. However, this barrier might be 
eliminated by the use of some pre-amplification procedure 
such as whole genome amplification (31). 

In summary, we have developed a universal labeling 
system that can significantly increase the multiplexing cap- 
ability for homogeneous detection. In addition, we have 
shown, using two examples, that the accuracy and robust- 
ness of the 2D label-based assays are comparable with 
those of well-established low-multiplex and low-throughput 
methods. Regarding the ease-of-use, low cost, inherent 
standardization and time and labor effectiveness of this 
approach, this system is expected to be used in screening 
pathogens responsible for many complex diseases and 
become a method of choice for high-throughput and multi- 
plex mutation/SNPs genotyping assays. Prospectively, the 
strategy will have the potential to upgrade current mid-plex 
assays from various complex solid-phase systems to a uni- 
versal homogeneous platform that can be performed on the 
widely available real-time PCR machines. 
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